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1.0  Introduction 

1.1  Background 

The  gust  front  defined  by  a sudden  wind  shift  and  temperature  change,  is 
a dramatic  consequence  of  the  thunderstorm  downdraft.  Such  a gust  front  is 
the  interface  between  the  warm,  moist  low-level  ambient  air  and  the  cool, 
nearly  saturated  downdraft  air  of  mid-level  origin.  The  horizontal  velocity 
of  the  mid-levels  contained  in  the  downdraft  momentum  is  additionally  accel- 
erated by  the  negative  (downward)  buoyancy  resulting  from  the  evaporation  of 
precipitation.  Tns  downdraft  is  deflected  into  a horizontal  plane  by  the 
ground  and  dissipated  by  interaction  with  surrounding  air  and  with  the  ground 
surface.  The  outflow  is  a potential  hazard  to  aircraft  due  to  wind  shear  and 
turbulence.  A previous  report  (Goff,  e^  ^. , 1977)  provides  initial  details. 
In  this  report  we  address  numerical  modeling,  mean  models,  turbulence,  flow 
patterns  and  auxiliary  observations. 

1.2  Objectives  and  Scope 

Several  objectives  are  primary  to  this  study.  These  are: 

(1)  investigate  and  categorize  the  turbulent  structure  of  the  subcloud 
layer  and  the  outflow  gust  front 

(2)  observe  aircraft  response  during  landing  approaches  through 
thunderstorm  gust  fronts  and 

(3)  investigate  suitable  numerical  models. 

' ' ! 

I Secondary  objectives  are: 

(a)  collection  of  wind  and  temperature  data  in  the  lower  levels 
during  gust  front  cases 

(b)  provide  a significant  data  bank  for  laboratory  and  computer 

i'y'it  I simulation  studies 


(c)  determine  quasi-steady  state  characteristics  of  thunderstorm 
outflow  and  gust  front  features  and 

(d)  investigate  use  of  Doppler  radar  and  satellites  for  gust  front 
detection. 

1 . 3 Research  Plan 


The  basic  organization  is  detailed  in  Federal  Aviation  Administration 
Report  No.  FAA-RD-77-119  and  is  not  here  repeated.  The  461  m tower  instru- 
mented at  seven  levels  to  444  m continued  in  operation  during  the  1977  pro- 
gram. Dual  Doppler  radars,  a WSR-57  weather  radar,  a surface  station  network 
and  an  F-4-C  aircraft  form  the  1977  research  observations'  base.  The  1977 
season  had  below  normal  thunderstorm  activity  and  no  exceptionally  strong 
gust  fronts  were  reported.  Significant  new  data  are  discussed  in  the  follow- 
ing sections. 

2.0  New  Data 


Several  interesting  gust  fronts  were  observed  jointly  through  the  use  of 
Doppler  radar,  aircraft,  tower,  and  satellite. 


(a)  (b) 

Figure  2-1  WSR-S?  weather  radar  display  for  (a)  1520:17  CST  and 
(b)  1550:15  CST  May  19^  1977,  0°  elevation.  Radar 

contours  at  approximately  10  dBZ  intervals.  Range 
marks  (circles)  are  at  40  km  intervals. 
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2.1  May  19.  1977 
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i 
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On  this  date  a large  squall  line  extending  more  than  300  km  (162  nm)  In 
a north-south  orientation  passed  through  Central  Oklahoma  (Figures  2-1a  and 
b);  Its  gust  front  embedded  In  light  precipitation  reached  the  tower  at 
1557  CST.  Figure  2-2  Is  a section  from  the  objective  analysis  of  the  tower 
data  covering  the  time  of  this  event.  Tower  data  were  filtered  before  plot- 
ting the  time-height  sections  (Goff,  1975). 


The  gust  frontal  zone  Is  characterized  by  moderate  shear  In  the  wind 
speed  component  normal  to  the  front.  A sharp  temperature  discontinuity  Is 
also  evident.  An  updraft  larger  than  4 m s"^  at  1557  CST  Is  followed  by  a 
downdraft  exceeding  2 m s~\  thus  creating  a somewhat  turbulent  zone  just 
behind  the  gust  front.  Little  surface  pressure  discontinuity  was  associated 
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Hgure  2~2  19  May  1977  KTVY-toioer  time-height  eeotione  of  etrecmlinest  poten- 

tial tempevature  and  wind.  Unite  are  m a~^  and  °K.  Streamlines 
are  drawn  using  wind  speed  normal  to  the  gust  front  (+  toward  the 
right) t wind  shear  parallel  to  the  front  (+  toward  low  pressure) 
and  vertical  velocity  (+  upward).  Thunderstorm  is  to  the  left; 
outflow  is  moving  from  left  to  right.  Each  analysis  is  450  m 
thick.  A time-to-epace  converted  1 km  distance  is  inidcated  at 
the  upper  right.  Rainfall  (rm)  trace  is  shown  at  bottom  of  figure. 
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Figure  2-3  Dual-Doppler 
wind  field  of  a gust 
front  19  May  1977^  for 
(a)  1526  CST  and  (b) 
1532  CST.  Wind  speeds 
are  indicated  by  length 
of  arrow.  Dashed  line 
indicates  gust  front. 
Tower  wind  at  444  m 
plotted  at  "+".  Col- 
trane  Rd.  site  surface 
wind  is  plotted  at  1 
and  Pennsylvania  Ave. 
site  surface  wind  is 
plotted  at  2.  Approach 
end  of  runway  at  Tinker 
Air  Force  Base  indi- 
cated by  ”A"  and 
approach  end  of  runway 
17R  at  Will  Rogers 
International  Airport 
by  "fl".  Full  barb  = 

5 m s-^;  half  barb  = 

2.5  m 8~^. 


with  the  gust  front.  Light  precipitation  began  a few  minutes  ahead  of  the 
front.  Figures  2-3a  and  b show  the  low  level  (center  of  beam  Is  250  m above 
ground)  dual -Doppler  derived  winds  at  1526  and  1532  CST.  Superimposed  are 
the  tower  winds  at  the  444  m level  and  the  surface  winds  at  the  Pennsylvania 
Avenue  (10  km  NW  of  tower)  and  Coltrane  Road  (9  km  SE  of  tower)  sites.  Note 
how  well  defined  Is  the  small  meso-low  at  the  north  end  of  the  gust  front  In 
Figure  2-3b.  This  characteristic  pattern  Is  similar  to  that  shown  by  Goff 
et^.  (1977). 

The  F-4-C  flying  at  460  m (1500  ft)  above  ground  level  (AGL)  just  north 
of  the  tower  on  east-west  flight  paths  measured  horizontal  and  vertical  winds 
and  turbulence. 

Horizontal  winds  (V)  are  calculated  using  the  true  airspeed  and  heading 

(plus  wind  vane  correction)  (V-_)  and  the  ground  speed  and  track  (V_,)  as 

ag  , . 9® 

determined  by  the  onboard  Inertial  navigation  system  (INS)  or 


Vertical  winds  (w)  are  computed  using 


" ■ ♦ “r/a  - V "Az"*  " -p.sC'  * L,® 

where  = true  air  speed 
a s angle  of  attack 
0 “ pitch  angle 
6^  * yaw  angle 

a^  « vertical  acceleration  of  aircraft 
'*n  “ vertical  motion  of  the  aircraft  at  time  T = 0 

e > pitch  rate 

Ly  * distance  from  accelerometer  to  angle  of  attack  measurement  point. 
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As  a measure  of  turbulence,  the  derived  gust  velocity  (U^^)  Is  used. 

This  value  Is  computed  from  aircraft  data  recorded  using  the  general  equation 


where  w 


weight  of  aircraft 

measured  Incremental  vertical  acceleration  for  normal 
equivalent  air  speed  at  sea  level 
gust  alleviation  factor 
air  density  at  sea  level 

change  In  lift  coefficient  with  angle  of  attack 
wing  area 


Corresponding  subjective  turbulence  evaluations  are  (as  given  by  NASA) 


light  Ujg  < 6.0  m s’^  (<  19  ft  s 

moderate  6.1  < ^ 9.1  m s"^  (20-29  ft  s“^) 

severe  9.2  s s 12.1  m s"^ (30-39  ft  s"^) 

extreme  U^g  > 12.2  m s"^  (^  40  ft  s"^) 


The  horizontal  winds  for  one  representative  pass  made  westbound  at 
1527:21-1528:40  CST  are  shown  In  Figure  2-4.  Only  light  turbulence  was 
reported  by  the  pilot  which  Is  In  agreement  with  the  recorded  turbulence 
and  vertical  winds  as  shown  In  Figure  2-5.  At  about  1530  CST  It  was  decided 
to  make  approaches  to  Tinker  Air  Field  (TIK)  (NE  of  Norman)  (Point  A In 
Figure  2-3a)  ahead  of  the  gust  front  and  then  to  make  an  approach  to  Will 
Rogers  International  Airport  (OKC)  (Point  B In  Figure  2-3a).  Figure  2-6  Is 
a plot  of  the  aircraft  recorded  winds  along  the  flight  path  during  an 
approach  to  TIK  at  1539-1541  CST.  The  winds  are  10  second  mean  winds  derived 
from  the  0.1  second  data  and  correspond  to  about  1 km  space  averaged  wind. 
Only  some  light  turbulence  and  no  precipitation  was  encountered.  The  air- 
craft was  then  directed  toward  an  approach  to  OKC.  The  gust  front  had 
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Figure  2-4  Westbound  winds  recorded  by  F-4-C  aircraft  1527:21-1528:40  CST 
May  19,  1977  at  460  m (1500  ft)  AGL. 
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Figure  2-5  Time  arose  section  of  F-4-C  recorded  derived  gust  velocities 

(A)  and  vertical  winds  (B)  for  1527:21-1528:40  CST  May  19,  1977. 

passed  OKC  about  1548  CST.  The  approach  was  started  at  1550  CST  and  the 
recorded  horizontal  winds  during  the  approach  are  shown  In  Figure  2-7.  Gusty 
lateral  winds,  turbulence,  lightning  and  some  light  precipitation  were  encoun- 
tered during  the  descent. 
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Figure  2-6  Winds  recorded  by  F-4-C  aircraft  during  approach  to  runway  at 


TIK  1539-1541  CST  May  19y  1977.  Numbers  in  parentheses  indi- 
cate altitude  above  ground. 


Figure  2-7  Winds  recorded  by  F-4-C  during  approach  to  runway  17R  at  OKC 
1550-1552  CST.  May  19,  1977.  Numbers  in  parentheses  indicate 
altitude  above  ground. 

Several  additional  approaches  were  made  to  TIK  ahead  of  the  gust  front. 
These  approaches  were  terminated  just  before  heavy  rain  crossed  the  airport 
and  as  fuel  became  critical.  The  approaches  encountered  conditions  similar 
to  the  first  one— no  turbulence,  no  precipitation.  This  gust  front  can  thus 


Figure  2-8  WSR-S7  weather  radar  scope 
presentation  for  1733  CST  June  12 y 
1977.  Radar  contours  at  approxi- 
mately 10  dBZ  intervals.  Range 
marks  (circles)  are  at  40  km  inter- 
vals. 


to 


Figure  2-9  Thunderstorm's  northeastern  edge  taken  from  the  ASD  F-4-C  on 
June  12y  1977  1733  CST  flying  at  1520  m (5000  ft)  MSL  showing 
sharp  leading  edge  of  storm  marking  separation  of  warm  air 
inflow  and  lower  cold  air  outflow. 
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F-Lgure  2-10  Same  atorm  as  Figure  2-9  except  central  portion. 

be  characterized  as  having  strong  low-level  southeasterly  wind  with  slight 
turbulence  prior  to  the  wind  shift  and  northerly  flow  behind  the  front.  This 
northerly  flow  is  turbulent  and  gusty.  Satellite  photographs  covering  this 
event  showed  a large  cirrus  cloud  shield  over  most  of  Oklahoma  which  obscured 
all  low  level  cloudiness  or  any  other  possible  identification  of  the  gust 
front  or  gust  frontal  activity. 

2.2  June  12.  1977 

An  extremely  unstable  air  mass  prevailed  over  Oklahoma  and  a low  pres- 
sure center  was  located  near  Altus,  Oklahoma  (ITS)  in  the  southwestern  part 
of  the  state.  Around  noon  (1228  and  at  1250  CST)  two  small  tornadoes  were 
reported  near  Cl inton-Sherman  AFB  [158  km  (85  nm)  west  of  NSSL].  These 
particular  storms  weakened  and  then  later  in  the  afternoon  formed,  along  with 


Figure  2-11  Same  storm  as  Figure  2-10  except  southern  portion. 


new  cells,  a north-south  squall  line  in  Central  Oklahoma  (Figure  2-8).  From 
near  1615  to  1800  CST  the  F-4-C  flew  a number  of  passes  just  above  the  gust 
front.  The  gust  front's  surface  position  was  visible  to  the  pilot  as  a line 
of  dust.  This  dust  was  not  sufficient  to  be  visible  in  the  photographs. 
Figure  2-9  to  2-11  are  representative  photographs.  Figure  2-12  shows  the 
turbulence  and  vertical  winds  encountered  while  Figure  2-13  has  a plot  of  the 
horizontal  winds.  The  corresponding  dual -Doppler  winds  are  also  shown  in 
Figure  2-13.  The  heavy  line  AB  is  the  aircraft  track  with  point  A corre- 
sponding to  the  time  1739:42  and  point  B the  aircraft  location  at  1744:30. 

The  line  A'B',  showing  the  aircraft  measured  horizontal  winds,  should  be 
superimposed  upon  AB  but  for  clarity  in  illustration  has  been  displaced  to 
the  right.  The  satellite  photo  (Figure  2-14)  taken  at  1600  CST  shows  the 
storm  with  some  indication  of  the  cold  air  outflow  (shown  by  arrow  in  Figure 
2-14). 


UDE  MAX  PLOT 
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Figure  2-12  Derived  gust  velocities  and  vevtical  wind  cross  sections  for 
flight  above  gust  front.  June  12,  1977  1739-1745  CST  at 
1520  m (5000  ft)  MSL. 
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Figure  2-13  Aircraft  observed  winds  and  dual-DoTpgler  winds  for  June  12, 
1977.  Contour  lines  are  at  10  dBZ  radar  reflectivity  factor 
intervals. 
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Figure  2-14  Satellite  photograph  taken  at  1600  CST  June  12,  1977.  Arrow 
indicates  entrance  of  cold  air  outflow. 

2.3  June  28.  1977 

On  June  28,  1977  strong  cold  air  outflow  associated  with  a squall  line 
reached  the  tower  at  1704  CST.  The  tower  analysis  covering  this  event  is 
presented  in  Appendix  A,  Case  I.  Surface  gusts  at  the  tower  reached  28.0 
m s'^  after  the  passage  of  the  gust  frontal  boundary.  A sharp  temperature 
discontinuity  is  evident  across  the  front  as  the  temperature  dropped  about 
6°C  in  5 min.  Pre-gust  front  updrafts  are  greater  than  6 m s~^. 

The  Norman  Doppler  real-time  display  of  the  velocity  and  reflectivity 
fields  at  1642  CST  are  shown  in  Figure  2-15.  At  this  time  the  leading  edge 
{ of  the  squall  line  was  about  10  km  away  from  the  tower.  The  reflectivity 
‘ pattern  (Figure  2-1 5a)  is  typical  with  weaker  values  along  the  edges  of  the 
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Figure  2-15  June  28^  1977  veal-time  Doppler  radar  display  of  (a)  radar 
reflectivity  pattern  with  reflectivity  factor  (dBZ)  scale 
given  at  right  and  (b)  Doppler  radial  velocity  with  velocity 
scale  (m  s"^)  at  right.  Range  marks  for  (a)  and  (b)  40  km. 
Elevation  angle  is  zero  degrees. 
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squall  line  and  with  numerous  embedded  cores.  The  central  core  values  are 
greater  than  40  dBZ.  The  velocity  display  (Figure  2-15b)  shows  clear  evi- 
dence of  strong  outflow  winds  ( > 32  m s”^  toward  the  radar)  along  the  forward 
edge  of  the  line  displaced  from  the  cores  by  an  appreciable  amount.  From  the 
reflectivity  display  alone  one  may  have  judged  the  center  portion  of  the 
squall  line  to  be  weaker  than  the  extremities.  The  Doppler  velocity  display, 
however,  shows  this  area  having  strong  winds  (gust  front).  This  Is  an 
example  of  how  Doppler  radar  can  detect  outflow  winds. 

Damage  from  this  squall  line  was  confined  to  electric  power  outages 
caused  by  wind  and  lightning  In  the  northeast  Oklahoma  City  area  about 
1730  CST.  Ihe  recording  system  at  the  tower  lost  power  at  1712  CST.  Between 
1930  and  2000  CST  hall  up  to  3.8  cm  In  diameter  fell  In  a rural  area  62  km 
south  of  the  tower. 

3.0  Mean  Gust  Front  Analysis 

To  examine  the  general  characteristics  of  the  gust  front  as  observed 
from  tower  measurements,  eight  of  the  1976  gust  front  cases  and  four  of  the 
1977  cases  were  averaged  to  give  a mean  analysis.  Only  gust  front  cases  when 
the  faster  recording  data  rate  (1.3  s and  1.7  s Intervals  In  1976  and  1.5  s 
Interval  In  1977)  was  In  operation  were  used  to  determine  this  mean  structure 
(see  Table  1). 

The  time-to-space  conversion, 

AX  = -cAt  (3.1) 

was  used  to  calculate  distance  relative  to  gust  frontal  passage.  Here  Ax  Is 
the  distance  traveled  by  the  frontal  boundary  during  a time  Interval  At  at  a 
speed  of  c.  A 50  x 10  array  was  determined  for  each  case  at  250  m Intervals 
In  the  horizontal  and  50  m Intervals  In  the  vertical.  Each  horizontal  grid 
point  represents  an  average  over  250  m.  The  grid  points  extended  from 
2.125  km  ahead  of  the  gust  front  to  10.125  km  behind  the  front  and  from  0 to 
450  m in  the  vertical  for  each  of  the  twelve  cases.  To  arrive  at  the 
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Table  1.  Oust  Frontal  Cases  used  in  Mean  Analysis 


Date 

Time  (CST) 

Speed  (m  s"^) 

Orientation 

(deg) 

10  May  1976 

0144 

14.2 

65 

12  May  1976 

0844 

6.6 

65 

22  May  1976 

2208 

16.3 

352 

26  May  1976 

0014 

13.8 

355 

26  May  1976 

0727 

13.3 

30 

29  May  1976 

2359 

5.0 

75 

30  May  1976 

1828 

6.1 

240 

23  June  1976 

2140 

13.0 

350 

5 May  1977 

1807 

6.0 

20 

16  May  1977 

2154 

7.8 

65 

19  May  1977 

1556 

9.8 

30 

26  May  1977 

2138 

5.6 

38 

* Orientation  refers  to  azimuthal  direction  of  the  storm's 
left  flank  with  respect  to  the  gust  front. 
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composite  diagram  (Figure  3-1)  the  twelve  cases  were  simply  averaged  point 
for  point. 

As  seen  In  Figure  3-1  the  mean  analysis,  extending  from  2.125  km  ahead 
of  the  front  to  5.125  km  behind  the  front.  Is  just  a smoothed  version  of  the 
Individual  cases.  Note  that  no  portion  of  the  zero  contour  of  the  relative 


wind  speed  component  Is  coincident  with  the  zero  on  the  distance  scale  as  one 
expects.  This  discrepancy  Is  probably  due  to  displacement  of  the  gust  front 
boundary  during  the  averaging  process  and  due  to  errors  In  the  speeds  of  the 
Individual  gust  fronts. 
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since  all  of  the  cases  used  In  this  composite  analysis  were  of  the 
squall  line-type  some  of  the  features  of  the  frontal  zone  for  the  squall  line 
as  discussed  by  Goff  ^ al^.  (1977)  are  present.  Namely,  the  updraft  In  the 
warm  air  ahead  of  the  front  Is  followed  by  a downdraft  behind  the  front. 

Here  the  downdrafts  are  considerably  weaker  (less  than  1 m s“^)  than  those 
observed  In  Individual  cases.  The  strength  of  the  downward  motions  were 
evidently  diminished  during  the  averaging  process  since  downdrafts  are  not 
consistently  observed  In  the  same  location  relative  to  the  front  as  Is  the 
case  of  the  warm  pre-gust  front  updrafts.  Likewise,  secondary  surges  are  not 
apparent. 

The  protruding  nose  feature  or  bulge  of  the  gust  front  Into  the  warm  air 
Is  elevated  about  100  m above  the  ground.  This  Is  In  agreement  with  Goff's 
et  (1977)  findings. 

4.0  Dynamical  Model  of  Gust  Front 

4.1  Gust  Front  Dynamics 


Previous  gust  front  analyses  have  been  limited  to  the  layer  below  the 
highest  level  of  the  Instruments  Installed  on  the  NSSL  tower  (about  450  m). 
Actual  gust  fronts  are  distinguishable  to  greater  heights. 


Numerical  simulation  may  be  used  to  partially  overcome  the  limitations 
of  data  availability  In  the  attempt  to  construct  an  accurate  description  of 
the  structure  of  the  thunderstorm  gust  front.  There  have  been  few  attempts 
to  simulate  actual  gust  fronts,  but  a similar  phenomenon  referred  to  In  fluid 
dynamics  literature  as  the  gravity  current  or  density  current  has  been 
studied  extensively. 
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The  work  by  Mitchell  and  Hovermale  (1977),  referred  to  hereafter  as  the 
M-H  simulation,  is  a simulation  of  a gust  front.  To  illustrate  the  evolution 
of  the  gust  front,  we  shall  discuss  this  model  in  some  detail.  They  used  a 
non-hydrostatic,  dry,  two-dimensional  primitive  equation  model  as  described 
by  the  following  equations: 


equations  of  motion 


thermodynamic  energy  equation 


t * “If  * *11  * “ p'h  " 


mass  continuity 


(4.4) 


where  x and  z are  the  horizontal  and  vertical  coordinates,  u and  w are  the 
velocity  components  of  the  x and  z coordinates,  p and  p are  the  pressure  and 
density  of  dry  air,  Cp  and  c^  are  specific  heat  capacities  of  dry  air  at 
constant  pressure  and  constant  volume. 


The  temperature  (T)  and  potential  temperature  (0)  are  obtained  from  the 
equation  of  state  for  dry  air  and  Poisson's  equation: 


and 


e = T(^)'^/^p  (4.6) 


where  R is  the  gas  constant  for  dry  air. 

The  terms  in  equations  (4.1)  through  (4.4)  represent  linear  diffusion 
expressed  for  a given  variable  y as  follows: 

2 2 

Fr(y)  = K H + K ^ (4.7) 


6 2 1 

where  and  K,  are  constant:  K = 2.25  x 10  cm  s"  and  K,  = 1.25  x 

6 2^  •!  i X Z 

10  cm^  s"  for  interior  points  and  = 0 at  the  boundaries. 


The  boundary  conditions  are:  (1)  for  w,  under  the  condition  of  flat 
bottom  and  top  surfaces. 


w 


(2)  for  u,  the  slip  condition  Is  used  at  the  top. 


s 0 at  z = H 
3z 


(4.9) 


and  at  the  bottom  surface  a "limited  slip"  formulation  Is  used,  that  Is, 


» 0 at  z = 0 
3z 


(4.10) 


and 


- •“Cn  I U^  I U^ 

Ijl)  = ° J — at  z = 175  m 

Az 


(4.11) 


where  t Is  the  stress  and  a non-dimensional  constant  that  represents  the 
average  drag  coefficient  over  the  lowest  175  m and  Uy  Is  u at  175  m. 

The  boundary  conditions  for  pressure  and  density  are  chosen  as 


|2.  « .pg  and  ^ (l^-)  = 0 at  z = 0,  H 


The  lateral  boundary  conditions  for  u,  w,  p and  p are. 


(4.12) 


u - 0.  U « 0.  If  = 0 and  If  « 0,  at  X = 0,  L.  (4.13) 


3x 


3x 


The  Initial  condition  used  for  the  experiment  Is  a horizontally  uniform 
Initial  vertical  temperature  profile,  ^^(z)  as  constructed  from  the  National 
Hall  Research  Experiment  sounding  data  which  were  obtained  near  Fort  Morgan, 
Colorado  on  a day  during  which  four  severe  convective  storms  occurred 
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locally.  The  Initial  pressure  Pq(z)  and  density  Pq(z)  are  computed  from 
Tq{z)  by  using  the  hydrostatic  equation  and  the  equation  of  state. 


and 


91nPo 

9z 


-=1 


!fo 

3x  " 3x  “ 3x 


(4.14) 

(4.15) 

(4.16) 


The  atmosphere  is  assumed  to  be  at  rest  initially. 
4.1.1  Downdraft  Simulation 


The  gust  front  results  from  a downdraft  developed  inside  a thunderstorm 
by  evaporation  and  the  drag  of  raindrops.  The  evaporation  cools  the  air  and 
increases  the  density  causing  a downward  acceleration.  The  drag  of  falling 
raindrops  also  produces  a downward  acceleration.  A complete  simulation  of  a 
realistic  downdraft  development  has  not  yet  been  done.  Mitchell  and  Hovermale 
assume  that  a cold  downdraft  somehow  develops  and  persists.  To  establish 
the  downdraft,  the  temperature  in  the  desired  downdraft  location  was  reduced 
each  time  step  by  a pre-determi ned  small  increment  during  the  first  half 
minute  of  the  prediction  period.  The  artificial  cooling  was  maximum  (approxi- 
mately 0.33*’K  per  time  step)  at  about  z « 1.4  km  and  x > 0.  After  the 
initial  cooling  during  the  first  half  minute,  the  potential  temperature 
field  is  shown  in  Figure  4-1.  Within  the  boxed  region  in  Figure  4-2a,  the 
potential  temperature  is  held  constant  for  the  remainder  of  the  forecast. 
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Figure  4-1  Potential  tempevature  field  (Q)  at  t = 0.40  minutes.  Contour 
values  range  from  306°  to  333°K  at  interval  of  1°K.  (Courtesy 
of  Mitchell  and  Hovermale^  1977.) 

4.1.2  Numerical  Computation 

The  governing  equations  (4.1)  through  (4.4)  are  numerically  integrated 
over  a 70  X 41  rectangular  grid  having  a horizontal  extent  of  16  km  ( x = 

233  m)  and  a vertical  extent  of  7 km  ( z = 175  m).  The  finite  difference 
formulation  of  equations  (4.1)  through  (4.4)  used  in  the  gust  front  simula- 
tion is  identical  to  that  of  Shuman  and  Hovermale  (1968).  In  order  to 
suppress  the  high-frequency  oscillations,  mostly  computational  noise,  a 
temporal  filter  described  by  Asselin  (1972)  was  applied  at  every  time  step. 
For  any  function  of  time  F(t),  the  filtered  function  F(t)  is  obtained  from 

- F(t)  + j [F(t-At)  - 2F(t)  + F(t+At)]  (4.17) 

where  v is  taken  as  0.5  in  this  simulation.  This  filter  has  a good  low-pass 
response  characteristic  and  does  an  effective  job  of  controlling  the  compu- 
tational mode  arising  from  the  centered- In- time  integration. 
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The  potential  temperature  field  at  t - 12.00  minutes  Is  shown  In 
Figure  4-2a.  The  contour  values  range  from  306“  to  333“K  at  the  interval  of 
TK  as  in  Figure  4-1.  It  Is  apparent  that  a cold  outflow  Is  formed  by  steady 
subsidence  of  cold  air  from  the  rectangular  domain  Inside  which  the  cold 
downdraft  is  artificially  maintained.  The  leading  edge  of  the  outflow  is 
characterized  by  a sharp  gradient  of  potential  temperature,  nearly  discon- 
ti nuous . 

For  the  purpose  of  illustration  only,  streamlines  are  depicted  as 
stream  function  contours.  The  stream  function  (j<(x,y)  is  defined  by 

u » - and  w = ||-  (4.18) 

and  calculated  by  solving  the  following  equation  by  a Liebmann  method 

0*0'  (i-  I?) 

The  flow  at  t = 12  minutes  is  depicted  by  Figure  4-3,  4-4  and  4-5.  The 

wind  blows  in  the  direction  of  the  streamlines  and  is  faster  where  stream- 
lines are  closer  together.  At  t = 12  minutes,  a strong  downdraft  of 

15  m s'^  Is  obtained  at  z = 1.4  to  2.0  km,  x = 0 within  the  Initially 

cooled  area  (Figure  4-1). 

The  downdraft  turns  Into  the  horizontal  as  It  approaches  the  ground, 
reaching  a maximum  horizontal  speed  of  24  m s~^  at  0.3  km,  x 3.5  and 
z • 0.7  km,  X » 8 km.  The  outflow  air  displaces  the  low-level  air  in  its 
path  toward  the  downwind  boundary.  As  this  air  approaches  the  boundary  it 
Is  artificially  forced  upward,  since  the  boundary  condition  allows  no  flow 
through  the  boundary. 

This  boundary  condition  and  the  Initial  rest  (no  motion)  condition  are 
significant  weaknesses  of  this  simulation.  Because  of  the  boundary  condi- 
tion, the  same  air  which  Is  forced  to  ascend  near  the  right  lateral  boundary 
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Figupe  4-2a  Potential  temperature  field  (Q)  at  t = 12,00  min.  Boxed  area 
denotes  region  of  fixed  0 values.  Contour  range  and  interval 
as  in  Figure  4.1  (Courtesy  of  Mitchell  and  Bovevmale,  19??). 
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must  become  part  of  the  downdraft  when  it  reaches  the  left  boundary.  However, 
the  integration  of  the  model  is  stopped  before  the  low-level  air  has  time  to 
reach  the  downdraft,  so  only  the  initially  high-level  air  enters  the  down- 
draft.  The  source  of  the  air  would  not  affect  its  ability  to  sustain  the 
downdraft  because  of  the  artificial  cooling.  The  undesirable  result  of  the 
boundary  condition  is  the  quasi-circular  circulation  pattern  depicted  in 
Figure  4-3a. 

A more  realistic  picture  would  include  flow  through  the  lateral  boun- 
daries and  an  initially  sheared  flow.  The  downdraft  would  then  be  fed  by 
mid-level  air  entering  the  domain  through  the  left  boundary.  If  the  low- 
level  air  had  initial  motion  toward  the  left  boundary,  a more  realistic 
updraft  tilted  in  an  up-shear  sense  might  result.  Of  course  cross-boundary 
flow  would  introduce  computational  difficulties.  Unfortunately,  in  this 
simulation  the  circulation  pattern  results  completely  from  the  introduction 
of  artificial  negative  buoyancy  in  an  artificially  closed  domain. 

Despite  these  weaknesses,  the  results  of  the  simulation  are  useful  if 
we  direct  attention  to  the  behavior  of  the  cold-air  surge.  By  comparing 
Figures  4-2a,  4-3a  and  4-4a  with  the  careful  analysis  of  Charba  (1974)  using 
data  from  the  NSSL  instrumented  tower.  Figures  4-2b,  4-3b  and  4-4b,  the 
similarity  is  very  striking.  Notice  the  difference  in  the  scales  of  the 
figures.  In  Figure  4-3  the  wavy  character  of  the  flow  within  the  cold  air 
is  evident  in  both  the  simulation  and  the  analysis,  especially  the  downward, 
then  upward  turn  immediately  upwind  from  the  gust  front. 

The  surface  wind  is  slower  than  the  horizontal  wind  at  the  higher 
levels  because  of  surface  drag.  The  drag  coefficient  (Cp)  used  for  the 
experiment  shown  in  Figures  4-1  through  4-5  is  0.02.  Mitchell  and  Hovermale 
ran  five  additional  cases  using  a different  value  of  Cjj  for  each  case.  Of 
course,  a stronger  surface  wind  occurs  with  smaller  values  of  C^.  Figure 
4-6  shows  the  time  variation  of  the  surface  wind  at  a fixed  point,  x » 

8.8  km.  The  surface  pressure  variation  and  the  time  of  the  cold  air  surge 
are  for  the  ^a;se  where  Cq  « 0.01.  It  is  interesting  to  note  that  the 
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Figure  4-Za  Streamfunotion  (^)  at  t = 12.00  min.  Contour  vnterval  is 

2.0  X 10^  8~^.  (In  all  streamfunotion  plots^  loibels  are 

scaled  by  -1.0  x 10~^  and  solid  contours  depict  aounter- 
alochoise  flow.)  (Courtesy  of  Mitchell  and  UoveTmalej  1977.) 
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Figure  4-3b  Same  as  Figure  4-2b  except  stream  function  analysis,  (Courtesy 
of  Charbat  1972.) 
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Figure  4-4a  Horizontal  velocity  field  (u)  at  t = 12.00  min.  The  evper- 
impoeed  dashed  (or  dotted)  line^  whe-yiever  it  is  given  in  any 
figure f represents  the  314°K  isentrope  in  the  corresponding 
0 field.  A and  B denote  surface  maxima  and  C denotes  the 
high  velocity  core  above  the  surface.  Contour  values  range 
from  -10  to  28  ms~^  at  intervals  of  2 ms~^.  (Solid  contours 
in  all  velocity  plots  represent  positive  values.)  (Courtesy 
of  Mitchell  and  Hovermale^  1977.) 
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Figure  4-4b  Same  as  Figure  4-2b  except  'Hnd  field  analyeie.  (Courtesy  of 
CharbOt  1972.) 
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Figure  4~S  Vertical  velocity  field  (w)  at  t = 12. 00  win,  A denotes  the 

relative  maximum  of  dcwmoard  motion  behind  the  front.  Contour 
values  range  from  -IS  to  10  ms~^  at  intervals  of  1 me~^. 
(Courtesy  of  Mitchell  and  Hovemalet  1977. ) 


Figure  4-6  Horizontal  velocity  (u)  as  a function  of  time  at  a fixed  point 
(x  ^ 8.8  kn)  at  the  surface  (z  0)  for  all  six  surface  drag 
cases.  For  the  Cp  = 1.01  oasst  the  arrow  denotes  the  onset  of 
cold  air  and  the  dashed  line  depicts  the  surface  pressure  vari- 
ation. The  amplitude  of  the  surface  pressure  rise  in  all  cases 
was  4,5  ± 0.3  ni>.  (Courtesy  of  Mitchell  and  Hovemalet  1977.) 
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Figure  4-7  Horizontal  velocity  (u)  as  a function  of  time  at  a fixed  point 
[ (x  9.7  km)  at  the  surface  (z  =>  0)  for  the  four  cases  given  in 

j Table  1.  Each  profile  is  labeled  with  the  maximum  surface  wind 

I speed  and  the  corresponding  maximum  potential  temperature  decrease 

j (LQ  at  z z SOO  m).  The  dashed  line  depicts  the  surface  0 vari- 

i ation  in  the  12**K  case.  (Courtesy  of  Mitchell  and  Hovermale^ 

I 1977.) 


Increases  of  the  horizontal  wind  speed  and  the  surface  pressure  started 
several  minutes  prior  to  the  cold  air  surge. 

Throughout  these  experiments  the  downdraft,  the  origin  of  the  gust 
front,  was  produced  by  an  artificial  reduction  of  temperature  Inside  the 
pre-determlned  area.  Of  course,  the  strength  of  the  resulting  downdraft  and 
gust  front  depends  strongly  on  the  amount  of  cooling  prescribed.  Figure 
4-7  shows  the  surface  winds  resulting  from  four  different  values  of  total 
cooling  at  500  m. 

Although  this  numerical  simulation  Is  quite  Incomplete  from  the  view- 
point of  the  thunderstorm  circulation.  It  has  provided  some  useful  quantita- 
tive Insights  Into  the  mechanics  of  thunderstorm  gust  fronts. 
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4.2  Origin  of  the  Gust  Front 


Our  knowledge  of  the  gust  front  structure  is  fairly  complete  after 
extensive  analyses  of  the  data  from  tower-mounted  meteorological  sensors, 
special  radar  and  radiosonde  observations,  and  aircraft.  The  thunderstorm 
outflow,  the  leading  edge  of  which  forms  a gust  front,  has  many  characteris- 
tics analogous  to  gravity  or  density  currents  from  both  laboratory  and 
theoretical  studies.  Me  have  just  seen  that  Mitchell  and  Hovermale's 
numerical  experiments  succeeded  in  simulating  many  aspects  of  the  gust 
fronts  by  using  artificial  cooling  to  generate  a downdraft.  This  mechanism 
of  gust  front  generation  is  essentially  the  Margules  Process  (1905). 

4.2.1  Marqules  Process 


Let  us  consider  as  an  initial  state,  two  air  columns  of  the  same 
height  (h)  standing  side  by  side  above  the  ground.  One  is  colder,  marked 
by  {!>,  and  the  other  is  warmer,  marked  by  {2}  in  Figure  4-8.  The  fraction 
of  the  area  {1}  to  the  total  area  is  3.  The  total  horizontal  dimension  is 
l and  hence  the  column  {1}  occupies  the  area  3£. 

The  pressure  of  the  air  at  the  top  of  each  column  is  p^.  The  tempera- 
tures at  the  top  of  {1}  and  {2}  are  Tj,^  and  Tj^2»  respectively.  Temperatures 
at  the  ground  are  denoted  by  Tq^  and  Tq2  and  the  pressure  by  Pg^  and  Pq2, 
respectively.  If  we  assume  that  both  columns  have  dry  adiabatic  lapse  rates, 
r^.  Tg^  and  Tg2  become 


and 


'01  'hi 


'02  'h2 


+ rjh 


(4.20) 


(4.20’) 


I 
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and  the  surface  pressures  are 


01 


’’h 

" ’hi 


and 


(4.21) 


p . p (IM) 


where  1 


= -E 


{4.2T) 


The  potential  energy,  PE,  of  a unit  cross  section  of  a column  is,  using 
the  hydrostatic  relation. 


PE  * / gpzdz  * 


p dz  - p^h 


(4.22) 


and  the  internal  energy  IE,  using  the  equation  of  state,  is 


C..  h , h 


IE  = C„  / pTdz  = X 0 P ^ 


Adding  (4.22)  and  (4.23),  the  total  energy,  E,  per  unit  cross  section 
of  a column  becomes 


E - PE  + IE  * / p dz  - p^  h. 

'■  0 


(4.24) 


Substituting  P ■ Pq  ( ’’y  ° -)  into  (4.24)  and  carrying  out  the  Integration, 
the  total  static  energy  ^r  unit  cross  section  becomes 
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Now  since  air  column  {1}1s  denser  than  column  {2},  It  may  be  presumed 
that,  lacking  any  restraining  mechanism,  column  {1}  will  seek  a position 
underneath  column  {2},  reducing  the  static  energy  and  resulting  In  the  con- 
figuration pictured  In  the  right  half  of  Figure  4-8.  If  we  compare  Equation 
4.26  to  a similar  equation  describing  the  total  static  energy  of  the  system 
after  the  rearrangement,  the  difference  will  represent  the  kinetic  energy 
generated  by  reducing  the  potential  energy  of  the  system. 

1'  p^h 
(4.27) 

The  unknown  temperatures  and  pressures  of  the  final  state,  as  well  as 
l'  and  p^  may  be  determined  from  the  conservation  of  the  masses  of  the  two 
columns.  Assuming  that 

P^i  “ Ph 
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and  therefore 


h'  = h , 


the  masses  of  the  two  columns  at  initial  and  final  times  are 
Integrating  the  hydrostatic  equation: 


^ g 


m^ 


; Pq-pJ. 


P02"Ph 


,.p:-Ph 


Since  m-j  - m^ 


sKPoi-Ph)  = ii'(p;-p;) 


and  because  m2  - m^ 


(1-6)  t (Po2-Ph)  ■ 6’  (P;-Ph)- 


Solving  (4.28)  and  (4.29)  fort'  and 


t* 


®^Poi“Ph^  * ^^■^^^Po2“Ph^ 


I 
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obtained  by 


(4.28) 


(4.29) 


(4.30) 


(4.31) 


p;  - Pi  - ^^pqi  “Ph><p6-Ph) 

^^Poi'Ph^  * ^^“®)^P02"Ph^ 

If  the  redistribution  of  the  air  columns  Is  simply  adiabatic,  then 


I ■ ^h2 

(4.32) 

i r -T 

I 'oi  ^hl  'pj,' 

(4.33) 

I 1 

Pi  ^ 

1 

(4.34) 

, 1 

T'  B T (—1^ 

ml  hi  'pj,^ 

(4.35) 

p6  • ®P01  * <’-®>P02 

(4.36) 

If  the  total  energy  of  the  system  Is  conserved  during  the  adiabatic 
redistribution,  and  If  the  system  was  Initially  at  rest  (KE  = 0),  then 

KE'  » TE  - TE'  (4.37) 


where  M Is  the  mass  of  the  system 


In  order  to  make  a sample  calculation  of  the  wind  speeds  which  might 
result  from  such  a process,  let  us  take  the  following  Initial  conditions: 

h - 5568.15  m 
Pq^  > 1000  mb 
Tq^  - 302.51“  K 
= 248.16*  K 

After  the  rearrangement. 


p^  > 500  mb 
Pq2  - 987.63  mb 
Tq2  = 307.51“  K 
Tj,2  * 253.16*  K 


p^  = 990.10  mb 
t'  ■ A 

p*  « 890.10  mb 
T^2  « 253.16*  K 
« 301.65*  K 
T|Jj2  “ 298.51®  K 
- 292.61*  K 

If  the  remaining  calculatlonstare  carried  out,  we  see  that  for  this 
hypothetical  situation,  which  we  might  view  as  a typical  case  of  evaporative 
cooling  In  a thunderstorm.  If  all  the  energy  was  converted  to  horizontal 
motion,  then  a typical  wind  speed  of  12  m s**^  could  result. 
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This  procedure  might  be  refined  to  produce  an  estimate  of  the  maximum 
surface  winds  to  be  expected  with  thunderstorms  If  the  thunderstorms  produced 
a gust  front.  The  temperature  of  the  colder  air  mass  could  be  taken  to  be 
the  wet-bulb  temperature  of  the  middle-level  air.  The  wanner  air  mass  could 
be  assigned  the  wet-bulb  temperature  of  the  low-level  air.  The  role  of 
vertical  momentum  transport  also  has  to  be  Included. 

4.2.2  Far-reaching  Gust  Fronts 

Gust  fronts  frequently  propagate  as  far  as  30-50  km  away  from  the  thunder 
storm  center  without  losing  much  of  their  kinetic  energy  (Figures  4-9  and 
4-10).  In  the  M-H  simulation  this  distance  was  typically  about  10  km  (Figure 
4-2a).  The  simulated  gust  front  leads  the  thunderstorm  by  a shorter  distance 
than  was  typical  In  the  case  study  of  Charba,  despite  the  fact  that  the 
maximum  horizontal  wind  speed  was  also  smaller  In  the  simulation  (28  m s~^ 
compared  to  32.5  m s~^  In  the  Charba  analysis),  and  also  despite  the  fact 
that  the  temperature  difference  between  the  downdraft  and  Its  environment  has 
about  the  same  magnitude  In  both  the  simulation  and  the  analysis  (compare 
Charba's  vertical  temperature  profiles  In  Figure  4-11  to  Figure  4-2a).  This 
weakness  of  the  simulation  has  not  been  adequately  explained. 

4.3  Microphysical  Processes  (One  Dimensional  Thunderstorm  Simulation) 

We  will  now  examine  the  microphysical  processes  which  cool  the  air  and 
create  thunderstorm  downdrafts.  The  cooling  Is  mainly  due  to  evaporation, 
and  the  downward  acceleration  of  the  air  Is  ge.^erated  by  the  negative  buoy- 
ancy of  the  colder  air  and  by  the  drag  of  falllnj  rain  drops. 

Figure  4-12  shows  the  scheme  of  cloud  microphysical  processes  used  by 
Ogura  and  Takahashi  (1971)  for  their  one  dimensional  (vertical)  simulation  of 
a thunderstorm.  The  microphysical  processes  considered  to  be  responsible  for 
changes  among  the  water  phases  (vapor,  liquid  and  Ice)  are  condensation, 
evaporation,  conversion,  sublimation,  glaciation  and  melting.  These  micro- 
physical processes  will  be  described  here  In  some  detail  and  their  effect  on 
the  gust  front  will  be  considered. 
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Figure  4-9 


Major  squall  line  features  at  2337  CST  May  31,  1969.  The  PPI 
radar  echo  trading  is  of  the  2x10  and  2x10^  nrfi  reflectivity 
factor  contours.  XY  is  the  section  line  for  the  plane  of  verti- 
cal cross  section  analysis  of  the  MKY  tMer  data  shotJn  in 
Figure  4-10.  AfBi,  A2B2t  etc.  are  section  lines  used  by  Charba 
(1968)  for  the  time-to-spaoe  converted  surface  network  data 
traces  of  the  stations  indicated.  ( Courtesy  of  Charba. ) 
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Fi^fUTB  4—10  A ocmpoBi/be  Bchetnatio  modBZ  combining  ifiB  feci'tux‘68  of  the  ccnolyzed  ond  deduced  etmotux^ 

of  the  windehift  and  gust  front  leading  the  aquall  line  of  May  31,  1969.  (Courtesy  Charba). 
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Figure  4-12  Scheme  of  cloud  nricrophyeical  prooeaeee  in  a thunderetom.  ( Courtesy  of  Ogura  aovi 

Takahoehi,  1971.) 


Ogura  and  Takahashi  parameterized  the  various  microphysical  processes  in 
an  attempt  to  evaluate  their  relative  importance.  The  condensation  rate. 


dt 


PI  =-^ 


dq 


vs 

dt 


for  0 >0 

My  - Hys 


(4.41) 


= 0 


for  q^<  q^3  ^ 


where  q^  and  q^^  are  the  mixing  ratio  of  water  vapor  and  its  saturated 
value,  respectively.  The  rate  of  conversion  from  cloud  droplets  to  rain 
drops , 

' ‘^’CONV  ' 


where  q and  q are  the  mixing  ratios  of  cloud  droplets  and  of  raindrops, 
c r 

respectively,  and  the  parameter  is  the  reciprocal  of  the  conversion  time, 
taking  values  in  the  range  of  0 to  0.2  s The  glaciation  rate 

= G q„  (4.43) 

” glac 


where  q^  is  the  mixing  ratio  of  ice  crystals  and  6 is  the  reciprocal  of  the 
glaciation  time  from  0 to  0.05  s"\  The  rate  of  sublimation  used  is 


P4  » 


7 X 10®  + 0.41  X 10^ 
«is 


(4.44) 


where  q^^  and  e^^  are  the  saturation  mixing  ratio  of  water  vapor  and  satura 
tion  vapor  pressure  measured  in  mb  over  a plane  ice  surface,  respectively. 
The  constant  f^j  is  used  in  the  computation  of  the  terminal  velocity  of  ice 
particles  and  takes  the  value  0.75  here.  The  rate  of  melting  is 


P5  ■ 2.27  X 10“®C  (T-273)(p^q^)°-®25p^-lf^-0.42  45j 
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C represents  the  ventilation  constant  given  by 


C » 1.6  + 0.57  X 10’3(V.)1-5  f -1 

1 0 


The  rate  of  evaporation  from  cloud  droplets  to  water  vapor 


-= 


do 

evap  evap 


(4.46) 


and  evaporation  was  assumed  to  be  instantaneous  as  long  as  the  air  was  not 
saturated  with  water  vapor.  The  rate  of  evaporation  of  raindrops 


P7  = 5.8  X 10-^(q  - q ) 
'^v  ^vs' 


(4.47) 


where  the  units  are  gm  gm’^  s"^  This  equation  is  a simplified  form  of 


, (^-1)  C(P3  q^)0*525 

P7  = - i — 

5.4  X 10°  . 


(4.48) 


where  e^^  is  the  saturation  vapor  pressure  over  a water  surface.  The  rate 
of  evaporation  of  ice  crystals  is 


‘’a  7 X 10®  + Ml  X 10_^ 


(4.49) 


which  is  similar  to  that  for  raindrops  except  for  different  values  of  the 
latent  heat  of  sublimation  and  vapor  pressure.  The  rate  of  evaporation  of 
melting  ice  crystals 
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P9  = - — — 


0.525  ^ -0.42 
0 


5.4  X 10®  + ^ 

®ws 


(4.50) 


These  parameterize tions  are  applied  within  the  framework  of  a one 
dimensional  model  describing  a right  circular  cylinder,  with  radial  symmetry 
and  coordinates  (r,  x,  z).  The  pressure  distribution  Is  not  considered  but 
Is  assumed  Identical  to  the  environmental  conditions.  The  equation  of 
vertical  motion  Is  simplified  to 
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(4.51) 


where  for  any  variable  A, 


A.  * ^ /^’^A  dx  at  r = a , 

a Zn  0 

T^  Is  the  virtual  temperature  and  the  added  subscript  zero  refers  to  envi- 
ronmental conditions.  The  first  term  of  the  right  hand  side  represents  the 
vertical  advectlon,  the  second  term  the  lateral  eddy  exchange,  the  third 
term  the  dynamical  entrainment  which  Is  required  to  satisfy  mass  continuity 
between  the  cloud  and  the  environment. 
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The  fourth  term  Is  the  buoyancy  and  the  last  term  the  drag  force,  assumed  to 
be  proportional  to  the  weight  of  cloud  droplets,  raindrops  and  Ice  crystals. 


Similarly,  the  thermodynamic  equation  Is  given  by 

2 

O'  a a'  a 


If  • <1?  ♦ '■d>  * 1*1  (vn  ♦ 7 “.(T-T.)  + 
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(4.54) 
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where  L„,  L.,  L.  are  the  latent  heats  of  evaporation  (600  cal  gm~  ),  subll- 

V S T 1 •t 

matlon  (680  cal  gm“  ) and  fusion  (80  cal  gm  ),  respectively.  The  continuity 
equations  for  water  vapor,  cloud  droplets,  raindrops  and  Ice  crystals  are 
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where  is  the  terminal  velocity  of  raindrops 
= 3.12  X 10^  cm  s“'’ 

and  is  the  terminal  velocity  of  ice  particles 


V,  = 3.12  X lo’  cm 


(4.59) 


(4.60) 


Ogura  and  Takahashi  conducted  an  initial  experiment  without  any  micro- 
physical processes.  Figure  4-13  shows  the  time-height  cross  sections  of 
vertical  velocity,  excess  temperature  and  liquid  water  content.  A strong 
updraft  forms  and  reaches  a steady  state  after  about  45  minutes.  No  down- 
draft  and  no  middle-level  temperature  deficit  occurs.  The  strong  tempera- 
ture deficit  appearing  at  11  km  is  apparently  caused  by  adiabatic  cooling 
in  the  upper  stable  layer  while  the  updraft  air  is  losing  its  momentum. 

A much  more  realistic  simulation  was  achieved  when  the  microphysics 
were  added,  as  shown  by  Figure  4-14,  which  is  similar  to  Figure  4-13  but 
the  concentrations  of  cloud  droplets,  raindrops  and  ice  crystals  are  also 
depicted.  The  maximum  downdraft  of  4 m s'^  developed  at  about  4 km  after 
about  65  min,  but  the  accompanying  temperature  is  hardly  significant.  The 
downdraft  starts  developing  first  in  the  lower  part  of  the  cloud  at  40  min 
and  spreads  to  the  higher  altitudes.  At  approximately  60  min,  the  down- 
draft  replaces  the  updraft  throughout  the  domain.  In  particular,  a strong 
downdraft  appears  at  the  melting  zone.  After  the  downdraft  reaches  the 
ground,  it  decreases  and  soon  dies  at  about  70  min.  The  maximum  updraft 
and  downdraft  throughout  the  life  cycle  is  17  m s"^  and  5 m s“\  respec- 
tively. 

The  excess  temperature  Increases  during  the  initial  40  min  in  a way 
similar  to  the  case  without  the  microphysical  processes.  After  that  time 
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*he  temperature  starts  decreasing  and  at  55  min,  a negative  temperature 
anomaly  developed  near  the  melting  level.  Figure  4-15  reveals  that  the 
maximum  downward  velocity  Is  observed  very  close  to  the  time  and  level 
where  four  processes  are  occurring:  evaporation  of  Ice  crystals  (P8)  and  of 
melting  of  Ice  crystals  (P9),  evaporation  of  raindrops  (P7),  and  melting  of 
Ice  crystals  (P5).  All  of  these  processes  appear  during  the  decaying  stage 
of  the  cumulus.  The  evaporation  processes  P7,  P8,  and  P9  have  the  same 
order  of  magnitude  while  the  melting  of  Ice  crystals  (P5)  cools  the  air  at 
a rate  about  one  order  larger.  The  total  cooling  Is  very  small  compared 
to  observations  in  real  thunderstorms,  largely  because  the  experiment  was 
one-dimensional,  precluding  the  possibility  of  considering  the  evaporation 
of  cloud  water  by  Invading  middle  level  dry  air. 

4.4  Dry  Air  Intrusion  (Two  Dimensional  Thunderstorm  Simulation 


Perhaps  the  weakest  feature  of  the  one  dimensional  simulation  Is  Its 
failure  to  portray  a strong  downdraft  with  accompanying  strong  evaporative 
cooling.  The  evaporation  process  was  hindered  by  the  absence  of  a source 
of  dry  air  into  which  water  could  easily  evaporate. 


Figure  4-16  Is  a schematic  representation  of  a typical  vertical  circu- 
lation pattern  within  a cumulonimbus  In  an  environment  with  vertical  shear 
(Ludlam,  1961).  Dry  middle  level  air  enters  the  storm  from  the  upwind  side 
and  Is  cooled  by  evaporating  precipitation  (Indicated  by  dashed  lines)  which 
falls  from  the  tilted  updraft.  The  resulting  Intense  cooling  causes  a 
downdraft  of  sufficient  magnitude  to  prevent  the  dry  air  from  destroying 
the  updraft. 

In  order  to  simulate  this  schematic,  we  need  at  least  a two  dimensional 
model.  Several  two  and  three  dimensional  experiments  with  vertical  shear 
have  been  performed  (Takeda,  1965,  1966,  1971;  Orville  and  Sloan,  1970; 

Hane,  1972;  Schlesinger,  1973,  1974;  Wllhelmson,  1974).  We  will  discuss 
Wane's  experiment  for  the  purpose  of  Illustration. 
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Figure  4-14  Time-height  oroae  eeotione  of  (A)  vertical  velocity,  (B)  exceee 
temperature,  (C)  liquid  and  eolid  water  content,  (D)  content  of 
cloud  droplets  (E)  content  of  rcdndropa  and  (F)  content  of  ice 
cryatale  for  a cloud  with  microphyaical  proceaaea;  Fp  ® 6.3®C/kin, 
Co  ■ 0.005.  (Courtesy  of  Ogura  and  Takdhaahi,  1971.) 
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T1ME(  MINUTES) 


TIME(  MINUTES) 


TIME)  MINUTES) 
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Figure  4-15  Time-height  aroaa  aeationa  of  the  ratea  of 
(A)  oondenaation  of  water  vapor  into  cloud 
drapleta  (PI)  and  evaporation  of  melting  ice 
oryatala  (P9);  (B)  evaporation  of  cloud  drop- 
leta  (P6)f  raindropa  (P7)  and  ice  oryatala 
(P8);  (C)  oonveraion  of  cloud  dropleta  into 
raindropa  (P2)  melting  of  ice  oryatala  (P5); 
and  (D)  glaciation  (P3)  and  growth  of  ice 
oryatala  by  aublimation  (P4)  in  unita  of 
10-7  g-1,  (Courteay  of  Ogura  and  Takahaahit  1971.) 
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Hane's  model  Includes  the  condensation  and  evaporation  processes,  liquid 
water  "In  bulk"  assuming  the  Marshall -Palmer  drop-size  distribution,  variable 
terminal  fall  speeds  of  rain  and  turbulent  mixing.  Not  Included  are  the  Ice 
phase,  horizontal  and  time  changes  of  air  density,  radlatlonal  effects  and 
the  air-ground  (sea)  Interaction. 


The  vorticity  equation  used  to  describe  the  flow  Is 

If  - -u  U -»  |a  - Ii(|i  t ||)-  -5^  + g ^ [||f|  (u"-u)] 

- W [ifl  (x”-"']  ' * ^ ['■x  U ] * ^ ['■z  If]  (4.61 ) 


where  the  vorticity  ® |j  " |j»  Is  the  virtual  temperature,  q^^  the  rain 
water  mixing  ratio,  q the  cloud  water  mixing  ratio,  C_  the  turbulent  mixing 
coefficient,  u"  and  w"  the  environmental  horizontal  and  vertical  wind  com- 
ponents which  vary  only  In  the  vertical  direction,  and  and  are  the 
time-space  dependent  horizontal  and  vertical  eddy  viscosity  coefficients. 


The  equation  of  continuity  Is 


^ (pu)  + --  (pw)  = 0 (4.62) 

where  P Is  air  density.  A stream  function  Is  defined  to  satisfy  the  continu 
Ity  equation. 


pu  *11  and  PW  = - 


The  vorticity  Is  expressed  by  the  stream  function  as 


pn 


1 3<|>  3p 
p 3Z  dz 


(4.63) 


(4.64) 
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The  first  law  of  thermodynamics  is  written 


where  T is  the  absolute  temperature,  the  dry  adiabatic  lapse  rate, 

the  latent  heating  or  cooling  rate,  T"  the  environmental  temperature, 
a ^function  of  z only,  and  ^9  the  constant  thermal  diffusivity  (500  m s ). 


The  prediction  equations  for  the  three  forms  of  moisture  are 


-I-  evaporation  - condensation 


(4.66) 


- autoconversion  - accretion  - evaporation  + condensation  (4.67) 
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+ autoconversion  + accretion  - evaporation  (4.68) 


In  these  equations  q^  is  the  water  vapor  mixing  ratio,  q^.  the  mixing  ratio 
of  cloud  droplets,  q^^  the  mixing  ratio  of  rain,  q"^  and  q"^  are  environmental 
quantities  that  vary  only  in  the  vertical  direction,  and  Vj  is  the  effective 
terminal  velocity  of  the  raindrop  distribution. 

The  autoconversion,  accretion,  evaporation,  condensation  and  terminal 
velocity  of  raindrops  are  evaluated  based  on  Kessler's  parameterization 
(Kessler  et  al_.  1961-64)  which  is  similar  to  that  used  by  Ogura  and 
Takahashi  (1971). 
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Hane  simulated  storm  development  for  various  vertical  wind  shear  cases. 

Figure  4-17  shows  the  result  of  the  moderate  shear  case  with  the  system  speed 
of  10  m s~^  at  t = 5.6  min,  which  Includes  the  two  dimensional  stream  func- 
tion,  (10  kg  m s ),  rain  water  mixing  ratio  (gm  kg”  ) and  cloud  water 
mixing  ratio  (solid  lines,  dashed  lines,  scalloped  lines  respectively)  In  the 
top  section,  temperature  anomaly  (°K)  In  the  middle  section  and  vertical 
wind  speed  (m  s~^)  In  the  lower  section.  The  negative  temperature  anomaly 
(In  the  middle  section  of  the  figure)  shows  a cold  dome-like  shape,  with  the 
coldest  air  at  the  ground.  The  flow  pattern  shows  that  the  downdraft  ori- 
ginated In  the  middle  levels  at  the  left-hand  side  of  the  storm,  comparing 
well  with  the  dry  Intrusion  depicted  In  Figure  4-16. 

The  downdraft  spreads  both  forward  (to  the  right)  and  to  the  rear  (left) 
near  the  ground.  There  Is  an  Inflow  from  the  right  to  left  near  the  ground 
that  sharply  slides  upward,  precipitating  and  producing  rain  (top  section  of 
Figure  4-17). 

The  outflow,  spreading  forward,  collides  with  the  low-level  Inflow 
underneath  the  central  portion  of  the  storm  (as  defined  by  the  cloud  water 
area  scalloped  In  Figure  4-17).  The  energy  dissipation  caused  by  the 
collision  may  retard  the  spreading  of  the  cold  air.  Hane  continued  the  time- 
integration  of  experiment  R1  up  to  t = 16.3  min.  The  outflow  did  not  move 
outside  the  storm  area  In  that  time.  In  experiment  R2  (strong  shear  with  the 
system  speed  of  19  m s"^),  the  leading  edge  of  the  outflow  at  t = 5 min  and 
15.7  min  never  leaves  the  storm  area  (Figure  4-18  for  t = 5 min  and  15.7  min). 

The  right  half  of  Figure  4-18  shows  that  the  typical  thunderstorm  circu 
latlon  pattern  seen  In  the  left  half  has  been  lost  by  15.7  min.  The  two 


dimensional  simulation  Is  a vast  Improvement  over  the  one  dimensional  model, 
but  Is  still  unable  to  capture  the  persistent  features  of  the  mature  thunder- 
storm with  a oust  front  propagatinq  30  to  50  km  ahead  of  the  storm. 
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Figure  4-16  Sdhematio  representation  of  the  motion  field 
within  a middle  latitude  aumuloninibuB  in  wind 
shear.  T indicates  height  of  tropopause.  The 
vertical  and  horizontal  scales  are  the  same 
(after  Ludlam,  1961). 

4.5  Twisted  Storm  CircuTation  (Three  Dimensional  Model) 

One  possible  reason  why  the  two  dimensional  simulations  fall  to  produce 
gust  fronts  which  propagate  far  away  from  the  storm  may  be  the  fact  that  the 
shear  between  the  low  level  flow  and  the  middle  level  flow  Is  exaggerated. 

In  two  dimensions  the  two  flows  are  exactly  opposing  each  other.  Typical 
thunderstorm  environments  have  somewhat  less  directional  shear  between  the 
low  and  middle  levels. 

It  may  well  be  that  a three  dimensional  simulation  with  a realistic 
shear  condition  would  produce  a gust  front  which  reaches  far  ahead  of  the 
storm.  In  recent  years,  some  effort  has  been  directed  toward  developing 
three  dimensional  thunderstorm  models,  with  limited  success  (Steiner,  1973; 

I Wllhelmson,  1974;  Schlesinger,  1975).  The  new  advanced  computers  should 

I expedite  these  efforts,  resulting  In  Increased  understanding  of  the  gust 

front  dynamics. 
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Figure 


12  16  2« 
DISTANCE  FROM  LEFT  BOUNDARY  (kin) 


Variable  fields  for  Experiment  R1  at  5.6  min.  Upper:  solid 
lines  are  streamfynation  [10^  kg  m~^  8~^)s  dashed  lines  rain- 
water mixing  ratio  (gm  kg~^);  scalloped  lines  oloudwater  mixing 
ratio  outline;  middle:  terrperature  anomaly  (°K);  lower:  verti- 
cal wind  speed  (m  s~^). 
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5.0  Sutiwary 


A large  data  bank  for  low  level  wind  shears  associated  with  gust  fronts 
has  been  compiled  from  observations  obtained  using  a 461  m tower  instrumented 
at  seven  levels.  These  data  are  on  magnetic  tape  and  are  available  for  use 
in  analytical  simulation  studies.  The  average  gust  front  is  characterized  by 
a 3 m s"^  updraft  immediately  above  the  gust  front  nose.  Mean  horizontal 
temperature  gradients  are  about  per  km  and  the  horizontal  wind  shear 
across  the  gust  front  averages  about  8 m s~^  per  km.  Gust  fronts  are  often 
characterized  by  secondary  surges  which  have  strong  shear  and  turbulence. 
Two-dimensional  models  have  been  developed  which  illustrate  the  thunderstorm 
downdraft  development  and  strengthening  due  to  negative  buoyancy  produced  by 
evaporation  of  raindrops  into  entrained  dry  middle-level  air.  These  down- 
drafts  which  can  exceed  8 m s“^  produce  a strong  outflow.  An  adequate  three- 
dimensional  model  still  needs  to  be  developed. 

The  gust  front  itself  is  a weather  hazard  because  of  windshear  and 
turbulence,  but  what  may  be  even  more  hazardous  are  the  windshears  and  down- 
drafts  directly  beneath  the  thunderstorm.  Here,  as  an  aircraft  passes  under 
the  thunderstorm,  the  aircraft  experiences  first  a strong  headwind,  then 
roughly  no  horizontal  wind  at  all  as  it  enters  the  downflow  area,  and 
finally  experiences  a strong  tail  wind.  These  wind  shifts  from  headwind  to 
tailwind  with  a vector  difference  of  40  m s"^  or  more  in  4 km  in  the  hori- 
zontal may  be  more  dangerous  than  the  well -organized  gust  front  and  also 
should  be  avoided  by  pilots. 

Detection  of  these  gust  fronts  has  been  a problem.  On  numerous  occa- 
sions it  has  been  all  too  obvious  that  surface  observation  stations  as  now 
distributed  form  a network  which  is  too  sparse  to  adequately  detect  (in 
sufficient  time  to  be  of  value  to  aviation)  all  gust  fronts  or  even  detect 
ones  having  strong  shear.  Supplemental  wind  observations  surrounding  an 
airport  offer  a viable  interim  method  for  thunderstorm  gust  front  location 
and  strength  determination  critical  to  aircraft  operation. 
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Remote  detection  by  satellite  is  not  reliable  due  to  (1)  cirrus  and 
other  high  and  middle  level  clouds  which  often  obscure  thunderstorm  bases  and 
that  volume  of  space  which  may  contain  the  gust  front;  (2)  inability  of 
weather  satellite  to  detect  "optically  clear"  gust  fronts  or  those  fronts 
which  contain  entrained  dust  particles  too  small  or  in  insufficient  numbers 
to  be  resolved  by  satellite  and  (3)  lack  of  resolution  in  the  infrared 
channels  used  during  the  night  time.  These  factors,  plus  the  limited  obser- 
vation frequency  now  available  or  planned,  severely  impact  satellite  usage  in 
gust  front  detection.  The  satellite  does  provide  another  means  besides  radar 
for  thunderstorm  area  determination. 

Since  a number  of  gust  fronts  are  in  optically  clear  air,  they  are  most 
often  not  observed  when  conventional  weather  radar  is  used.  The  radar  will 
indicate  precipitation  intensity  but  gives  no  indication  of  the  presence  or 
absence  of  a gust  front.  Doppler  radar  [Goff,  et  ^. , 1977]  can  det.ct  these 
gust  fronts  and  provides  the  most  promising  observation  tool. 

6.0  Conclusion  and  Recommendation 

Gust  fronts  and  associated  thunderstorms  are  hazardous  to  aircraft  opera- 
tion. These  fronts  are  extremely  variable  as  to  location  and  strength  and  are 
difficult  to  forecast  even  though  a general  life  cycle  is  present  consisting  of 
(1)  initial  downdraft  turning  into  an  outward  moving  cold  air  mass  closely  cen- 
tered near  the  downdraft  to  (2)  an  interface  possibly  more  than  ten  miles  in 
advance  of  the  storm  before  dissipating.  While  low  clouds  and  dust  occa- 
sionally mark  the  leading  edge  of  the  outflow,  usually  the  gust  front  is 
optically  invisible  and  thus  rarely  detected  by  conventional  weather  radar  or 
satellite  observation. 

Due  to  the  character  of  gust  fronts,  surface  observations  will  detect 
the  phenomena  and  adequate  spacing  of  wind  observations  in  the  vicinity  of  an 
airport  can  provide  an  interim  solution.  For  longer  range  detection  required 
in  air  traffic  control,  microwave  Doppler  radar  offers  the  highest  potential 
and  this  option  should  be  actively  pursued. 
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Development  of  an  acceptable  forecast  procedure  appears  to  be  some  time 
off.  While  simplified  techniques  can  provide  an  estimate  of  the  maximum 
Intensity  based  on  environmental  winds  and  temperature,  they  are  not  adequate 
for  determining  which  thunderstorms  have  gust  fronts  nor  where  gust  fronts 
are  in  relationship  to  the  thunderstorm.  Therefore  observation  and  tracking 
of  the  gust  front  appears  to  be  the  best  short  range  "forecast"  option. 

It  1s  recommended  that  continued  emphasis  and  major  support  be  placed  on 
development  of  both  ground  and  airborne  sensors  which  will  alert  the  pilot 
and  contr'oller  to  wind  shears.  Modeling  efforts  should  continue  with  the 
realization  that  It  may  be  some  time  before  an  adequate  forecasting  procedure 
Is  developed. 


APPENDIX  A 

Gust  Front  Cases,  1977 

CASE 

DATE 

TIME  (CST 

A 

5 May 

1807 

B 

13  May 

1346 

C 

16  May 

2154 

D 

19  May 

1556 

E 

26  May 

2138 

F 

28  May 

1035 

G 

31  May 

0746 

H 

12  Jun 

1513 

I 

28  Jun 

1704 

Cold  Front  Cases 

J 

10  Dec  1976 

1021 

K 

28  Jan  1977 

0311 

FIGURE  LEGEND 


Left  Pages: 

Time-height  sections  of  streamline  analysis  and  vertical  velocity 
(m  s~^).  potential  temperature  (°K),  horizontal  wind  speed  component  parallel 
to  the  gust  front  (m  s~^)  and  wind  speed  component  normal  and  relative  to  the 
gust  front  (m  s~^).  Each  objective  analysis  is  450  m thick;  time  increases 
from  right  to  left.  Outflow  is  moving  from  left  to  right.  Time-to-space 
converted  1 km  distance  is  indicated  in  the  upper  right.  Rainfall  and  pres- 
sure change  traces  are  shown  at  the  bottom  of  figures  when  necessary. 

Right  Pages: 

WSR-57  surveillance  radar  display  with  echo  contouring  at  0°  elevation. 
Radar  contours  are  at  approximately  10  dBZ  intervals.  Range  marks  (circles) 
are  at  40  km  intervals.  Time  clock  (CST)  is  in  the  upper  right.  The  KTVY-TV 
tower  is  located  at  358°  and  38  km. 

Remarks  consist  of  information  taken  from  Storm  Data  (U.S.  Department  of 
Conmerce)  and  personal  observations. 


Figure  A-A-1 


Figure  A-A-2 

Case  A: 

Date:  5 May  77 

Time  of  gust  front:  1807  CST 

Speed:  6.0  m s"^ 

Orientation:  20“ 

Pressure  jump:  no  pressure  rise 
Rain  began:  1850  CST  ended:  2239  CST 
Total  rainfall:  8.9  mm 

Maximum  intensity:  23.9  mm/hr  1921-1929  CST 

Remarks:  Gust  front  accompanying  a large  squall  line.  Although  gust  frontal 
winds  were  weak  in  the  tower  vicinity,  wind  gusts  up  to  35  m s“l  had  been 
reported  earlier  near  Mangum,  Oklahoma  (about  200  km  WSW  of  the  tower).  At 
1931  CST  an  apparent  small  tornado  touched  down  in  the  northeast  part  of 
Duncan  nearly  130  km  SSW  of  tower. 

The  wind  shear  associated  with  the  component  of  the  wind  normal  to  the  front 
is  weak  as  well  as  the  temperature  discontinuity  across  the  front.  Of  inter- 
est, though,  is  the  strong  updraft  (greater  than  8 m s"'  in  the  top  tower 
layer)  prior  to  frontal  passage. 
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Figure  A-B-2 
Case  B: 

Date:  13  May  77 

Time  of  gust  front:  1346  CST 

Speed:  8.3  m s"^ 

Orientation:  340° 

Pressure  jump:  no  pressure  rise 

Rain  began:  1430  CST  ended:  1451  CST 

Total  rainfall:  1.1  mm 

Maximum  intensity:  2.7  mm/hr  1430-1451  CST 

Remarks:  Weak  outflow  accompanying  a small  line  of  thunderstorms.  Wind  shear 
and  temperature  discontinuity  are  not  strong  across  the  frontal  boundary  but 
an  updraft  in  excess  of  4 m s"’  occurs  just  ahead  of  the  front.  The  strongest 
shear  in  the  wind  component  normal  to  the  front  began  at  1402  CST.  A second 
drop  in  temperature  as  well  as  an  increase  in  upward  vertical  velocity 
occurred  at  this  time  showing  evidence  of  a secondary  surge.  No  severe 
weather  was  associated  with  this  line  of  storms. 
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Figure  A-C-1. 


Figure  A-C-2 


Case  C: 


Date:  16  May  77 

Time  of  gust  front:  2154  CST 

Speed:  7.8  m s"^ 

Orientation:  65° 

Pressure  jump:  1.8  mb  2115-2159  CST 
Rain  began:  2347  CST  ended:  0130  17  May 
Total  rainfall:  14.5  mm 

Maximum  intensity:  24.8  mm/hr  2347-0017  17  May 

Remarks:  Moderately  strong  outflow  case  associated  with  a dissipating  squall 
line.  The  gust  front  was  identified  on  the  WSR-57  radar  by  a "thin  line." 
The  strongest  horizontal  shear  and  temperature  discontinuity  are  within  the 
gust  front  region.  A turbulent  zone  extends  from  just  ahead  of  the  front  to 


m s 


2157  CST  with  updrafts  greater  than  4 

3ms"'.  The  pressure  analysis  shows  a gradual  pressure  rise  of  1 


and  downdrafts  in  excess  of 

mb  in 

about  10  minutes  beginning  at  2146  CST.  Measurable  rainfall  did  not  begin 
at  the  tower  until  2347  CST  about  two  hours  after  the  wind  shift. 

This  line  had  produced  tornadoes  earlier  in  extreme  western  Oklahoma  but  no 
severe  weather  in  the  vicinity  of  the  tower  was  reported. 
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Figure  A-D-2 
Case  D: 

Date:  19  May  77 

Time  of  gust  front:  1556  CST 

Speed:  9.8  m s“^ 

Orientation:  30° 

Pressure  jump:  no  pressure  rise 

Rain  began:  1551  CST  ended:  2150  CST 

Total  rainfall:  52.0  mm 

Maximum  intensity:  unknown,  some  data  missing  due  to  power  failure 

Remarks:  Moderately  strong  outflow  associated  with  a large  squall  line.  The 
most  intense  portion  of  the  line  remained  south  of  the  tower.  The  related 
gust  front  which  was  imbedded  in  light  precipitation  reached  the  tower  at 
1556  CST.  The  gust  frontal  zone  is  characterized  by  moderate  shear  in  the 
wind  speed  component  normal  to  the  front.  A sharp  temperature  drop  is  also 
evident.  An  updraft  above  4 m s“'  at  1557  CST  is  followed  by  a downdraft 
nrpat.pr  than  ? m c-1 . thus  creating  a somewhat  turbulent  zone  along  the 


greater  than  2ms 
frontal  boundary. 


A short  narrow  tornado  occurred  north  of  Foster  (112  km  S of  tower)  at 
1630  CST.  One  tornado  took  place  at  1723  CST  near  White  Bead  (94  km  SSE  of 
tower)  and  one  reported  close  to  Paoli  (88  km  SSE  of  tower)  at  1736. 
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Case  E: 

Date:  26  May  77 

Time  of  gust  front:  2138  CST 

Speed:  5.6  m s"^ 

Orientation:  38° 

Pressure  jump:  no  pressure  rise 

Rain  began:  2309  26  May  ended:  0128  27  May 

Total  rainfall:  28.9  mm 

Maximum  intensity:  52.2  mm/hr  2309-2332  CST 

Remarks:  Gust  front  associated  with  a large  squall  line.  There  is  moderate 
wind  shear  in  the  wind  component  normal  to  the  front  but  a somewhat  weak 
temperature  discontinuity  near  the  surface.  Downdrafts  are  persistent  from 
the  2148  to  2153  CST.  Downdrafts  of  about  4 m s"’  extend  down  to  200  m above 
the  ground  at  2149  CST. 

Wind  speeds  of  35  m s"l  were  measured  near  Mustang,  Oklahoma,  at  2135  CST. 
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Figure  A-P-1, 


Figure  A-F-2 
Case  F: 

Date:  28  May  77 

Time  of  gust  front:  1035  CST 

Speed:  8.7  m s“^ 

Orientation:  155° 

Pressure  jump:  .3  mb  1033-1036  CST 
Total  Rainfall:  None 

Remarks:  The  satellite  photo  at  1000  CST  depicted  a thunderstorm  induced 
mesoscale  cold  front  as  a well-defined  arc  of  cumulus  clouds  from  south  of 
Stillwater  to  about  50  miles  west  of  Enid.  Its  movement  was  southward  at 
approximately  10  m s~'. 

At  the  tower  wind  speeds  gusted  to  13  m s"^  with  the  arrival  of  this  boundary 
and  the  wind  direction  shifted  from  a southerly  to  east-northeasterly  direc- 
tion. No  strong  wird  shear  or  turbulence  was  associated  with  this  boundary 
and  there  was  no  severe  weather  within  the  tower  vicinity.  The  pressure  rose 
si ightly. 
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Figure  A-G-l. 


Figure  A-G-2 


Case  G: 

Date:  31  May  77 

Time  of  gust  front:  0746  CST 

Speed:  6.5  m s"^ 

Orientation:  85° 

Pressure  jump:  1.3  mb  0744-0750  CST 
Total  Rainfall:  None 

Remarks:  Boundary  accompanying  a roll  cloud  reached  the  tower  at  0746  CST. 
There  is  moderate  horizontal  wind  shear  and  turbulence  along  and  behind  this 
boundary.  An  updraft  in  excess  of  3.0  m s"l  is  just  ahead  of  the  gust  front 
followed  by  a downdraft  greater  than  2 m s"^ . 

^ A significant  pressure  change  was  associated  with  this  disturbance.  A pres- 

sure increase  of  1.0  mb  took  place  in  3 minutes  beginning  at  0744.  No 
f measurable  precipitation  accompanied  this  disturbance. 


Figure  A-H-2 
Case  H: 

Date:  12  Jun  77 

Time  of  gust  front:  1513  CST 

Speed:  4.2  m s"^ 

Orientation:  60° 

Pressure  jump:  .2  mb  1511-1514  CST 
Total  Rainfall:  None 

Remarks:  Weak  outflow  case  associated  with  small  line  of  thunderstorms. 

Weak  wind  shear  and  temperature  discontinuity  are  evident  in  the  gust  frontal 
zone.  The  most  intense  part  of  the  line  was  south  of  the  tower. 
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Figure 


(WSR-57  photo  unavailable) 


Figure  A-I-2 
Case  I: 

Date:  28  Jun  77 

Time  of  gust  front:  1704  CST 

Speed:  7.8  m s"^ 

Orientation:  47® 

Pressure  jump:  no  pressure  rise 
Rain  began:  1711  CST  ended:  1808  CST 
Total  rainfall:  5.1  mm 

Maximum  Intensity:  unknown,  data  loss  due  to  power  failure 
(See  additional  remarks  in  text.) 


(WSR-57  photo  unavailable) 


Case  J: 

Date:  10  Dec  76 

Time  of  cold  front:  1021  CST 

Speed:  8.0  m s"^ 

Orientation:  60° 

Pressure  jump:  .5  mb  1020-1022  CST 

Remarks:  Strong  cold  front  which  passed  the  tower  at  1021  CST.  The  front 
was  accompanied  by  a sharp  temperature  discontinuity.  The  temperature 
dropped  8°C  in  4.minutes.  Also  accompanying  the  front  was  an  updraft  in 
excess  of  4 m s"'  and  moderate  shear  in  the  horizontal  component  normal  to 
the  front.  Pressure  analysis  shows  a .5  mb  rise  across  the  frontal  surface. 
The  analysis  shows  characteristics  very  similar  to  that  of  the  gust  front 
analysis. 
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Figure  A-K-1. 


(WSR>57  photo  unavailable) 


Case  K: 

Date:  28  Jan  77 

Time  of  cold  front:  0311  CST 

Speed:  12.0  m s"^ 

Orientation:  80° 

Pressure  jump:  .5  mb  0310-0313  CST 

Remarks:  Cold  frontal  passage  at  the  tower  marked  by  a slight  temperature 
fall  at  the  surface  followed  by  a gradual  temperature  Increase.  A strong 
Inversion  In  the  tower  layer  Is  evident.  A couplet  of  up  and  down  motions 
existed  across  the  frontal  boundary  with  updrafts  In  excess  of  4 m s"'  ahead 
of  the  boundary  and  downdrafts  greater  than  4 m s“l  behind  It.  The  wind 
speeds  Increased  from  nearly  calm  up  to  15  m s~'  at  the  surface.  A gradual 
pressure  rise  of  .5  mb  was  associated  with  frontal  passage  as  In  the  10  Dec  76 
case. 


84 


APPENDIX  B 


TABLE  DESCRIPTION 


Each  table  describes  the  number  of  tower  observations  made  per  hour  per 
day  for  one  month.  Time  (Central  Standard  Time)  Increases  from  left  to  right 
and  dates  increase  downward.  For  example,  in  Table  B1  for  the  hour  beginning 
at  1800  CST  on  the  16th  day  of  the  month,  362  observations  were  made  and 
recorded  on  magnetic  tape. 
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= XX88XXXXR88XR88X8lRRXRXR88l8X8X 
R X88XX8XX  8XRXXX8XlX8XXRXXBlX88X 
8 R8X8XX88iX8XXXX8RXX888RX8«iXX8X 
^ 8 X X X X X 8 X X X 8 X R 8 R R X 8 X X X R 8 8 X X S I X 8 X 8 

8 8 8 X R 8 X X R X X 8 8 8 X X 8 8 8 X 8 X R 8 8 8 1X888 

8 X888X8XX88888X88I8888X888  188X8 

8 flIX  RXXXX8XXX8XXIX8R88XX8  IXXXR 
8 XlX  X X X X X X 8 X 8 X 8 8 I X X 8 X R 8 X 8 1 8 8 8 8 
8 X I X S 8 8 8 8 X 8 8 8 8 X X 8 I 8 X 8 8 8 X 8 8 IR888 
8 X I 8 8 X X X X X 8 8 8 X 8 X 8 I 8 8 8 8 X 8 8 8 1 8 8 X8 
8 X I X 8 8 X X X X X 8 X 8 X 8 8 I 8 8 X 8 8 X X X 1 8 X 8 8 
8 RlRX88X8XXX8XXR8lXX88X8888l8888 
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IXKISIMIIIIXKMMiiVIII 
SiSXISliKIIIIMIXIMIIKI 
KfllliXIKSIKXKIilliliiK 
KXXSSillKKXiXXKIiSIIII 
XXXSXNXXXXKlIRKXXIRXXi 
RRRRRRRRRRRRRRRRiRRRR 
RRRRRRRRRRRRRiRRRRiRR 
RRRRRRRRRRRRRRRRRRRRR 
RRRRRRRRRRRRRRRRRiRRR 
RRRRRRRRRRRRRRRRRRRii 
RRRRRSRRRRRRRRiRRiRiR 
RRRiRRRRRiRRRilRRRRRR 
RRRRRRRRiRRERRiRRRRRR 
RRRRRRR  RRRR  RRRRsRRR 
RRRRRR  RRRRRSRRiRRiR 
RRRSilRisRRRRRRRRRRRRR 
RRKRRXRRRRRRRRRiRiRRR 
RRRRRRRRRRRRRRRRRRRR 
RRRRRRRRRRRRRRRRRRRRR 
RRRRRRRRRRRRRRRRRRRRR 
RRRRRRRRRRRRRRRRRRRRR 
RRRRRRRRRRRRRRRRRRRRR 
RRRRRRRRRRRRRRRRRRRRR 
RRRRRRRRRRRRRRRRRRRRR 


R R R R R 
R R R R R 
R R R R R 
R R R R R 
R R R R R 
R i R i R 
R R R i R 
R R R R R 
R R R R R 
R B R R R 
R R « R 
R R « R 
R R R I 
R R R R 
R R R R R 
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IXIIIIIIIIIIIIIIIiSlilllilllSXB 
SllllllllllllllilllllllllllBllli 
I X I I 8 I I i I I I I I I I 8 I 8 I X X I X B B B B B X B 
X X X X X X X X X X X X X 8 X 8 X X X B X X X X X B B X B X 
X X X X X X X X X X X X X 8 X X X X X B X X X B 8 8 B X B X 
XXXBXX8XXX8XBBXcBXX8XBXXBXX8X8X 
XXBXX8BXXBKXX8XXXXBXXBXXXX8BBXB 
BBXX8BBBBBBXXBXXXXXBB8XBXXBXX8B 
BX8XX8XXXXX8BXXBXBXBBXXXXX8BBXB 
XXBXX8BBBBBXBXX8BB8BXBXXXX8BBBX 
XBBXBXBBXXBXBBBX2BSKBBXXX8XBX8B 
XX.BSBXBBXBsBBXXBBmXXBBXBXXXXB- 
XBB8X8BBXXK&BXX8XXBBX8XB88BXXXX 
X 8 B 8 X X X X 8 B X B B B B X R 8 B B B s B XX8BBBB 
B B B S 8 8 B X X B B X B B X X X B X X t;  8 8 X X B X X B B 
BBBBXXBBBXXBXXBBXBBKXBXXXXXXXBX 
XXBBBBXX8BBiBKXXBBBXXBXXXXXXXKX 
BBX88XXBBBXBBBXXBXBBBBX88XBXXBB 
KBX8BXBBBXBBXBBXBBBBKXXBBBXXBXB 
XXBXBB8BBB8XBBBBXBXBBXBXBBBXXXB 
BBBBXXBBBXBXBBBXXBBXXXXIXXXBXBX 
XXXBB8BBBBXBXBX8BBXBXBXXXBXXXXX 
XBBBXBXXXBXXXXBXXXXXXXXXBXXXXXX 
BXBXXXBBXXXBXX8BXXXXXXXXXXXXXXK 
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s 8X8  8 X 8 8 8 8 8 8 8 X 8 8 X X 8 88888X88X8 

£ 8XX  8 8 8 X 8 8 X 8 8 8 8 8 8 X X 88X8888888 
a X 8 8 8 X 8 X 8 X 8 a 8 8 8 X 8 X X 8 8 8 8 8 8 8 X 8 8 8 

a 8 X X 8 8 8 X 8 a X 8 8 X 8 8 X X X 8 R8XRX8X888 

a 8 8 8 8 X 8 X 8 8 8 8 8 X 8 8 8 8 X X 8888888X88 
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